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We show that the electron-phonon coupling in graphene, in contrast with the non-relativistic two- 
dimensional electron gas, leads to shifts in the phonon frequencies that are non-trivial functions of 
the electronic density. These shifts can be measured directly in Raman spectroscopy. We show that 
depending whether the chemical potential is smaller (larger) than half of the phonon frequency, the 
frequency shift can negative (positive) relative to the neutral case (when the chemical potential is 
at the Dirac point), respectively. We show that the use of the static response function to calculate 
these shifts is incorrect and leads always to phonon softening. In samples with many layers, we find 
a shift proportional to the carrier concentration, and a splitting of the phonon frequencies if the 
charge is not homogeneously distributed. We also discuss the effects of edges in the problem. 

PACS numbers: 78.30.Na, 81.05.Uw, 63.20.Kr 



I. INTRODUCTION 

The discovery of graphene, a thermodynamically stable 
two-dimensional (2D) crystal, whose electronic proper- 
ties, described in terms of a half-filled 7r-clectronic band 
with Dirac electrons, can be controlled externally, has 
stirred great interest in the scientific community since 
the demonstration of a theoretically predicted 2 ^ anoma- 
lous integer quantum Hall effect 4 ' 5 . Unlike other 2D 
electronic systems, such as MOSFET heterostructures, 
graphene is easily accessible to optical probes. Further- 
more, in contrast to ordinary semiconductors where the 
different types of disorder can be distinguished through 
the temperature dependence of the transport properties^, 
graphene does not show any strong temperature or mag- 
netic field dependence in its electronic transport 7 - that 
allows an easy discrimination between different types of 
impurities. Hence, local probes such as scanning tun- 
neling microscopy (STM) and single electron transistor 
probes, will play a fundamental role in the understand- 
ing of nature the effects of disorder in graphene-based 
systems. 

Raman spectroscopy has been one of the most 
successful experimental methods used to study these 
systems ^ 10 ' 11 . In particular, it has been shown that it 
is possible to measure the number of graphene layers on a 
Si02 substrate with great accuracy, leading to an efficient 
and fast method to characterize graphene in situ. An- 
other interesting feature of these measurements is that, 
even for a single graphene layer, the phonon frequency 
measured in Raman shifts by a few wavenumbers, from 
point to point in space£&i£. Moreover, the observation 
of a D-line, which is Raman forbidden in translational 
invariant graphene, indicates the presence of disorder in 
the samples at electronic scale. 

We show that this Raman shift can be associated 
with the earlier experimental evidence for charge inho- 
mogeneity in undoped, unbiased, graphen o 7 ! 12 . There- 



fore, Raman spectroscopy can be used to map the dis- 
order in graphene layers, and hence, help to shed light 
on the nature of the disorder scattering in these materi- 
als. The understanding of the nature of impurity scat- 
tering in graphene is fundamental not only for the devel- 
opment of electronic devices based on carbon, but also 
may help to solve theoretical puzzles such as the discrep- 
ancy found between the theoretically predicted universal 
value of the conductivity 1 —, 4e 2 /(nh), and its experimen- 
tally observed^ value of 4e 2 //i (the so-called "mystery 
of the missing 7r" ) , and the absence of weak- localization 
effects 7 - (a topic that has generated intense theoretical 
debatei^^^^^^). 

In this paper we show that the shift in the phonon 
frequency in graphene has its origin on the polarization 
of the electrons due to the ion motion. Since graphene 
is a perfect hybrid between a metal and a semiconduc- 
tor there are two contributions to the polarization func- 
tion: one comes from intra-band transitions and another 
that originates on inter-band transitions. We show that 
the simplest approximation based on the static response 
is incorrect and predicts a reduction of the phonon fre- 
quency (softening of the lattice). The correct dynamic 
response is used to calculate the phonon frequency shift 
and it is shown that the phonon frequency can either de- 
crease (softening) or increase (hardening) depending on 
whether the phonon frequency is either larger or smaller 
than twice the chemical potential, respectively. We also 
show that the intra-band dynamic response vanishes at 
long wavelengths in a translationally invariant graphene 
sheet, while the inter-band contribution is finite. Never- 
theless, in disordered graphene we expect the intra-band 
contribution to be of the order of the inter-band one, in- 
dicating that disorder is important for the measurement 
of the Raman shift in graphene. 

The paper is organized as follows: in section |TT] we 
present the model for the electrons, phonons, and their 
coupling in graphene; section IIIII discusses the problem 
of the shift of the phonon frequency due to the electronic 
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polarization in graphcne and we consider both the static 
and the dynamic response; in section ITVl we examine the 
problem of phonon frequency shifts in bilayers and mul- 
tilayers within the same framework; section fVl contains a 
discussion of the problem of edges in finite samples and 
also the main conclusions of our work. We have also 
included one appendix with the details of an analytical 
model for the in-plane phonon modes in graphene, and 
also discuss the effect of defects and edges in the phonon 
spectra. 



II. THE MODEL. 

In the absence of disorder the Hamiltonian for electrons 
and phonons in graphene can be written as H = He + 
Hp +He~p, where (we use units such that h = 1 = ks), 

He = -t J2 ( r V- + h - c -) ~ M J2 c U Ca >* > M 

i,a=A,B 

is the free electron Hamiltonian, where [i is the chemical 
potential, c a> i (c' a 4 ) annihilates (creates) and electron on 



sublattice a = A, B on site Ri in the honeycomb lattice 
(spin indices are omitted throughout the paper), and to « 
2.7 eV, is the nearest neighbor hopping energy. Hp is the 
phonon Hamiltonian: 



Hp = J>, 



b 



(2) 



where b^ a {b\ a ) annihilates (creates) a phonon with mo- 
mentum q, and i — TA , LA , TO , LO, are the four 
phonon modest. In the following, we focus on the trans- 
verse optical (TO) modes near the T and K and K' points 
of the Brillouin zone. The TO band shows little disper- 
sion with a frequency ojo ~ 0.19 eV. 

We assume that the electron-phonon coupling arises 
from the modulation by the phonons of the Carbon- 
Carbon distance, a — 1.42 A, which leads to a change 
in the nearest neighbor hopping to- The dependence of 
to on distance / has been extensively studie d 22 ' 23 : 



dto/dl 



a ' 



6.4 eVA" 



(3) 



The resulting electron-phonon interaction is: 



He-p = (dh/dl) ckcsk+q {.TA q [l - e*( k+ ^72 - e l ( k +^ b /2] - x BcL [l 

k,q 

- (V3/2)y Bq [e ik ' a - e ik ' b ] } + h. 



72- 



ikb 



/2] 



(V3/2)y M [e 



i(k+q)-a _ *(k+q)-b 



(4) 
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where a and b are the unit vectors of the honeycomb lat- 
tice, and x a q, y a q (a = A, B) are given by the polarization 
of the phonon of wavevector q. They can be written as: 




y/2M C 0}q 




(5) 



where Mq — 1.2 x 10 4 m e is the carbon mass (m e is the 
electron mass), and the vector (ot\, a2, «3, CK4) is normal- 
ized to one. 

In order to obtain the polarizability of the TO mode, 
we use a central force model (see appendix [3} which 
leads to a phonon dispersion which can be calculated 
analytically 24 . This model is adapted from similar mod- 
els for tetrahedrally bonded lattices 2 ^. The details of the 
model are described in the appendix \^ where it is il- 
lustrated by some applications. The honeycomb lattice, 
even in the limit when the bonds are incompressible, can 
have shear deformations, leading to a vanishing shear 
modulus. Because of it, the model shows a flat transverse 
acoustical branch at zero energy. The optical modes, on 
the other hand, induce significant changes in the bond 
lengths. We focus here on a single optical mode, whose 



energy we take from experiments. The polarization is 
fixed by symmetry considerations. Hence, the model is 
needed only to describe the coupling to the electrons. 
The only coupling consistent with nearest neighbor tight 
binding model used to describe the it bands is the one 
that we are using. 

The polarization of the non-degenerate mode at the K 
point in the Brillouin zone is: 

(ai, a 2 , a 3 , a 4 ) - (1/2, i/2, -1/2, i/2) , (6) 

and we have a doubly degenerate mode (a Dirac phonon) 
with polarizations: 



(ai, 0-2,013,0:4) = (l/V2,-i/V2,0,0\ , 

(o x , 02,03,04) = (0, 0, 1/V2, 1/V2) • (7) 

For comparison, the polarization of the two optical modes 
at the r point can be written as: 

(01,02,03,04) = fl/v^O, -1/V2,0^ , 
(01,02,03,04) = (0,1/^,0,-1/^) • (8) 
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K+k+q 




K'+k 



FIG. 1: Diagram which describes the modification of the 
phonon propagator due to electron-hole transitions. See text 
for details. 



III. SINGLE LAYER GRAPHENE 

We are interested in the modification induced by elec- 
tronic transitions of the frequency of a phonon with 
wavevector Q. The electronic transitions which describe 
these processes are given, approximately, by: 



7Yq = 3a/2y^ c^ k cbk' 

k 



k (XAQ - X B Q + 



iVAQ + iVbq) + h.c. . 



(9) 



where we assume that the main contribution arises from 
transitions close to the Fermi level. In this limit, we can 
use the continuum limit and expand the energy of the 
electrons around the K and K' points, leading to the 
Hamiltonian: 
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where v-p = Stoa/2 w 6 eV A is the Fermi-Dirac velocity, 
and </>k = a,rctaxi(k y / k x ) is the angle in momentum space. 

A typical diagram which describes the renormalization 
of the phonon propagator in second order perturbation 
theory is given in Fig. [I] . The convolution of electronic 
Green's functions shown in the diagram is formally iden- 
tical to the charge susceptibility of grapheme: 



d 2 kdij' [G° A c Z(k, w')G^ pty (k + q, w 



G^(k,c')G c B T y (k + q^ + c')] 



(11) 



where G°£ c (k, u) 



(cr pty (k^)) 



with 



A, B is 



the electronic Green's function for the occupied (empty) 
states. The main difference between (| 1 1 1) and the charge 
susceptibility of graphene is that the charge susceptibility 
includes an overlap factor which suppresses completely 
transitions between the valence and conduction band at 
q = 0. 

From Fig. Q] we can immediately obtain the shift in the 
phonon frequency due to the polarization of the graphene 
layer due to particle-hole excitations: 



5u> 



Q 



27V3V 



ai ~ a 2 + 1013 + ta^l 



dl 



x(uQ,q^ 0) , 
(12) 



where we have used that Q — 3\/3a 2 /2 is the area of the 
unit cell. 

Graphene, from the electronic point of view, is a hybrid 
between a metal and a semiconductor: the polarization 
involves not only inter-band excitations (as in the case of 
the ordinary electron gas) but also intra-band excitations 
(as in the case of a semiconductor). The full susceptibil- 
ity x can be separated into an intra- and an inter-band 
contributions: 



X (^)= X intcI X9,w)+X intra (<7^)- 



(13) 



The intra-band contribution to the susceptibility was 
originally calculated by K.W.-K. Shung for graphene in 
re f . [2(1 and more recently it has appeared on refs. [27l . 



A. Static approximation 

A commonly used approximation on the electron gas 
problem is to replace the dynamical response x( tJ — 
uoQ,q — 0) by the static one: x( w = 0><? = 0). This 
approximation is usually justified in ordinary metals be- 
cause the Fermi energy fi is much larger than the phonon 
frequency so that the phonons respond to a time averaged 
electron distribution. In graphene, however, this is not 
necessarily so. In what follows we will study the effect 
of an static response and compare it with what happens 
when a dynamic response is calculated instead. We will 
show that these two approximations give very different 
results. 

At finite doping, the compressibility sum rule for the 
charge susceptibility leads to the equation: 



Iim n X intra (q,w = 0) = -D(n) 



2/i 



(14) 



where D{^l) is the density of states at the Fermi level. 
The number of carriers per unit cell, uq , is given by, 



3\/3 
2tt 



(fcpa) 



(15) 



The coupling to optical phonons involves terms which 
couple the two sublattices A and B, eq.©, while 
the charge susceptibility leading to eq. (fi"4"l) is the re- 
sponse to the total charge density operator: p(q) = 

Sk( c !ik+q c ^k + c sk+ q c sk). This change modifies the 
overlap factors, leading to an angular factor, when q — > 0, 
equal to sin 2 (6>) or cos 2 (0), depending on the polarization 
of the phonon. The average of this term gives rise to a 
factor 1/2 with respect to the (dia gonal) charge suscep- 
tibility calculated in refs. mH^m^H^ ' 
In addition, we have: 



<(q ^0,^ = 0) = - 



(16) 
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where the angular factor, as for the intraband suscep- 
tibility, depends on the polarization of the phonon and 
whose average is always 1/2. 

The main contribution to the integral in eq. (|16[) comes 
from k ~ A, where A is the high energy cut-off, so that 
this expression depends on details of the bands at high 
energies away from the Dirac point. Nevertheless, the 
change of the susceptibility with electronic density is in- 
dependent of the cut-off and can be readily calculated: 

lim (x intc » - x inter (0)) = 

* x k*«( g ->0,« = 0)«--V (17) 

which is of the same magnitude as the intra-band shift, 
eq.([T4]). 

Inserting JUJ, {T7P and fTS]) into (HI]), we find, in ad- 
dition to a density independent shift: 



Sluo = 7= -^r 



dto 
dl 



,V2 



McVQto 



(18) 



Expressing Sluq (in eV), and replacing uq by the density 
per unit area, n (expressed in cm -2 ), we find: 



Sto Q (eV) w -3 x 10- 9 n 1/2 (cnr 2 ) , 



(19) 



is the expression for the shift of the phonon frequency in 
the static approximation. For typical electron (or hole) 
densities, n w 10 11 — 10 12 cm~ 2 , the shifts are of the 
order of a few wavelengths (or degrees Kelvin), within 
experimental accuracy. 

Let us first notice that this result indicates that there is 
a decrease of the phonon frequency, that is, a softening of 
the lattice. This result is generically expected on physical 
grounds since a high density of electrons leads to the 
screening of the ion-ion interactions, reducing the elastic 
coupling in the lattice, and hence leading to a softening 
of the phonons. 



B. Dynamic approximation 

The real part of the intra-band susceptibility is given 
by^: 



5R[ X lntr >,g^0)] 



2irui u) 



2/i 1 



In 



2/i — LU 



2/i + lu 



(20) 



This expression is rather different from the static result 
(fT4|) . We note that the limits of lu = and q — ► with 
lu — > and q = do not commute. Moreover, we clearly 
see that (f20| changes behavior whether w/(2/i) is smaller 
or larger than one, and the susceptibility has a logarith- 
mic singularity in lu = 2/i. 

For 2^ > w a self consistent calculation of the po- 
larizability shows the existence of the two dimensional 
plasmon, which needs to be taken into account. We find: 



K[ X mtr ><2/i,g^0)] 



-q 2 H 



(«)] 



(21) 



where Lu p i(q) — (2e 2 /iq)/eo is the plasmon frequency 
(e is the electric charge and eo the dielectric constant of 
graphene). For 2/i <C lu we find: 



3?[ X intr >» 2/^^Q)] 



9 1 



3tt 



(22) 



Observe the change of sign in the expression of the sus- 
ceptibility in the two limits. More importantly, one can 
clearly see that these expressions vanish when q — > 0. 
This effect occurs because at q = the states associated 
with these transitions are orthogonal. Hence, in a system 
with translational invariance the intra-band transitions 
give no contribution. Nevertheless, in the presence of dis- 
order (or a finite sample), the electron mean free path, £, 
(or the system size, L) acts naively as a infrared cut-off 
and one would expect to see a non-zero effect. Replacing 
([2"Tj) and (|2"2")l into (TT2"j) , and assuming that wq ^> tu p i, we 
find: 



6lu^(2^>lu q ) 



£u;jf ra (2/i < loq) 



_fdto\ 



n 



1/2 



vpq 

UJQ 



2V6^V dl J MctOQto 

( tpvpg 



V2 \dl 



ii 



AI, 



C^Qto I 



(23) 



(24) 



Q 



Expressing <5wq in eV, and replacing uq by the density 
per unit area, n, expressed in cm -2 , we find: 

+ 1.4 x 10' 6 n 1/2 (qa) 2 , (25) 

-1.0 x l(T 18 n 3/2 (ga) 2 . (26) 



^ tra (2/i > u Q ) 
fc^ tra (2 M < WQ ) 



We stress, once again, that this shift vanishes as q — » 
and hence there should be no shift in the phonon fre- 
quency in a translationally invariant graphene sheet. 
Nevertheless, in the presence of disorder this is not nec- 
essarily the case. 

In order to include disorder in the calculation one 
would have to dress the fermion propagators in Fig. [T] 
by disorder and include vertex corrections to that dia- 
gram. These calculations are beyond the scope of this 
paper. Instead, we will follow a naive approach and sim- 
ply introduce a cut-off in q of the order of the inverse of 
the electron mean free path, £, which is known to be of 
order of 0.1/xm in these systems 7 . For typical electron 
(or hole) densities, n w 10 12 — 10 13 cm -2 , and I ~ 10 3 a 
\l « 0.1/im) the shifts are of the order of 10 -6 - 10~ 5 eV. 
For these concentrations and wavelengths, the plasmon 
frequency is co p i w 0.01 — 0.04eV, so that the assumptions 
leading to eq. are justified. Notice that while for large 
doping, /i > luq/2, the intra-band contribution leads to a 
hardening of the phonon, for low doping, /i < loq/2 there 
is softening of the phonon mode which depends directly 
on the amount of disorder in the system. 

The inter-band susceptibility is: 



1 



^[X inter (9^0,w)]« V / kdk 



4ej 



"lu 2 -4e 2 



(27) 



As in the static case, we find a large contribution which 
is independent of the carrier concentration and depend 
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on the high-energy cut-off. As before, we consider only 
the density dependent contribution: 



<( ff -0, W )]«— ^-^log 



lo — 2/i 



to + 2/i 



.(28) 



The first term in this expression reproduces the static 
limit of ^ lntor . The second term gives a correction which 
is more important for lo ~ 2/i, and cancels the static 
contribution as wq/fi — * oo. Comparing (|20[) with (|28p 
we find that: 

3?[ X intra (c, g = w/vf)] = 5W ntcr K 9 = 0)], (29) 



implying, from (|26p . that: 



jointer 



q (2/i > wq) 

int 

Q 



<W£ ter (2/z < 



+3 x 10~V /2 , (30) 
-2 x l(T 21 n 3/2 . (31) 



Once again, for /i > loq/2, the inter-band contribution 
leads to a hardening of the phonon, while for /i < loq/2 
there is softening of the phonon mode. Notice, that the 
numerical value of the inter-band contribution is small 
for densities of order 10 12 cm~ 2 when compared with 
the intra-band contribution estimated in the presence of 
disorder. It may well be that in disordered graphcnc 
the intra-band transitions dominate over the inter-band 
transitions. Hence, the final result may vary with the 
amount of disorder in the samples. 



IV. BILAYERS 



The same suppression applies to inter-band transitions, 
as, in any case, the modulation of the hopping involves 
transitions from the A sublattice to the B sublattice. 
Using dimensional arguments, which are also valid for 
single layer graphene, at the neutrality point we find, 

Xbulyerfa' w = °) K I^PA-L > and > at nnite fillings, 



lim XbUaycr(q^ = 0) OC fc 2 /^ 
|q|-»0 



(33) 



As in the single layer case, the inter-band susceptibility 
includes a contribution determined by the high energy 
cut-off, and density dependent term, which, also on di- 
mensional grounds, depends on density as the intra-band 
susceptibility, (|33| . 

For the inter-band contribution, using the reduction in 
the amplitude at the A sublattice mentioned earlier, we 
find: 



X&r(q,w = 0) 



t 2 , 



(34) 



where Abu = tx/vp is a high momentum cut-off above 
which the assumption that v-pk -C t± ceases to be valid. 
The integral in this expression has a logarithmic de- 
pendence on Abe, similar to the logarithmic divergences 
which characterize the charge susceptibility of a bilayer—. 

As in the case of the single layer, when we insert eg . (j3~4"|) 
into the expression for the shift in the phonon frequency, 
we find an term which is independent of the number of 
carriers, given by k-p. Taking it out, and neglecting log- 
arithmic corrections, we find: 



The previous analysis can be extended to a bilayer sys- 
tem. For simplicity, we consider here the static limit 
only. The model for the in-plane phonons considered 
before needs no changes. The shift in the phonon fre- 
quency is also given by the electronic susceptibility shown 
in Fig. PQ, and given in eci. tfTTj) . 

In a bilayer, however, the wavefunctions corresponding 
to the low energy electronic states have a small ampli- 
tude, aAk ~ k|k|/tj_ in the orbitals hybridized through 
the hopping t± with an orbital in the next layer (t± « 0.3 
eV is the inter-layer hopping energy). The relevant sus- 
ceptibility involves a convolution of the Green's function 
of sites in both sublattices, so that the reduction in low 
energy spectral weight at the sites with a neighbor in the 
next layer will reduce the susceptibility. 

The amplitude of an electronic wavefunction at energy 
^plkl 2 /^ <C 1 is of order a# ~ V~ x / 2 , where V 
is the volume of the system, on sites of the sublattice not 
connected to the second layer, defined as sublattice B. 
The amplitude on sites in sublattice A, where the sites are 
connected to the second layer is of order a a ~ UF|k|/tj_ x 
as- Then, the contribution to the susceptibility from the 
low energy electron-hole pairs is of the form: 



c bilayer 

ou>q oc 



Xb ^ er (q,u, = 0)« / dk 



\k(k + q) n k+q - ?i k 
t± e k+ q — £k 



dl J McuJQtx 



(35) 



(32) 



We expect a similar dependence for other multilayer sys- 
tems, as the main ingredient in this estimate, the changes 
in the low energy density of states in the two sublattices 
in each graphene layer, is independent of the number of 
layers in the stack. If the carrier density differs signif- 
icantly among the layers^ 2 -, we expect that phonons at 
each layer will be shifted by a quantity which depends 
on the local charge. Note that a crossover to a shift typ- 
ical of single layer graphcnc will take place at loq ~ t±. 
A similar crossover will occur if fj, ~ t±. 

A bilayer system can show a gap in the electronic spec- 
trum, when an applied field or chemical doping breaks the 
symmetry between the two layers. In this case, the elec- 
tronic states close to the gap are mainly localized in one 
of the layers. The polarizability shown in the diagram 
in Fig. [T] acquires a layer index, and is different in the 
two layers. Hence, we expect that the in-plane phonons 
in each layer experience a different frequency shift. In 
a first approximation, the phonons in the layer where 
the states at the Fermi energy have highest weight show 
the largest shift. Using dimensional arguments similar to 
those leading to eq. (f3"2")) . we expect that, when the Fermi 
wavevector is much smaller than the wavevector at the 
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center of the band of the biased bilayer, kp <C &o ~ A/up, 
the factor n ~ fcp in eci. ([3"5)) is replaced by ^ oc A 2 . 

The full dynamical response of a bilayer under a per- 
pendicular applied field is quite complex when the chem- 
ical potential is close to the gap edges^, with an anoma- 
lously large imaginary part. Hence, low energy phonons 
in a biased bilayer should be significantly damped. 



V. CONCLUSIONS 

We have analyzed the effect of a finite concentration 
of carriers on the frequency shift of phonons in electri- 
cally doped graphene samples. We have not considered 
changes due to modifications of the force constants, as- 
sociated to distortions of the a bonds. The analysis pre- 
sented here shows that the shift in the optical phonon 
frequencies in electrically doped graphene samples can 
be observed in Raman experiments^ ' 10 ' 11 , and it can be 
used to estimate the carrier density, or alternatively, the 
strength of the electron-phonon coupling. Notice that in 
an ordinary 2D electron gas the density of states (and 
charge susceptibility) depends only on the electronic ef- 
fective mass m* and is independent of the electronic den- 
sity. Therefore, for an ordinary 2D electron gas the fre- 
quency shift is essentially uniform and independent of 
disorder. For Dirac fermions, however, because of the ef- 
fective Lorentz invariance in the continuum limit, we can 
write an equivalent of Einstein's relation between energy 
and mass: /i = m*v^ (where the Fermi-Dirac velocity 
now plays the role of speed of light), indicating that the 
effective mass is energy dependent and vanishes at the 
Dirac point (/i = 0). Therefore, the effect described here 
does not work in an ordinary 2D electron gas. 

Another interesting consequence of equation (|24[) is 
that the Raman shift should be larger close to extended 
effects such as edges, dislocations and cracks 8,9,10 . The 
reason for that is the so-called self-doping effect discussed 
in great detail in ref. @: because of the poor screening 
properties of Dirac fermions, the Coulomb interactions 
remain long-ranged and an electrostatic potential builds 
up at the edges of the system, shifting the position of 
the surface states, and reducing the charge transfer to or 
from them. In this case the system, in order to maintain 
charge neutrality, can transfer charge to/from extended 
defects. This charging transfer is only halted when the 
charging energy of the edges is compensated by the ki- 
netic energy of the electrons. Thus, extra charge and a 
large density of states can be found at the edges of sam- 
ples. In this case, according to eq. (|24[) . the Raman shift 
should change as a function of the distance from the sam- 
ple edges (being larger at the edge). We have estimated 
that for edges of size 0.1-1 /i m the charge transfer is or- 
der of 10 -4 -10 -5 electrons per carbon (5n « 10 11 — 10 12 
electrons per cm -2 ) and hence the Raman shift is also of 
the order of a few wavenumbers but slightly larger than 
the effect produced by bulk disorder. We also notice that 
this effect is not possible in the ordinary 2D electron gas 



because screening leads to a uniform charge distribution. 

The shift of phonons with energies comparable or 
larger than the Fermi energy is determined by the dy- 
namic electronic response function, which is significantly 
different from the static one. In this regime, the shift 
changes when /i rs luq/2, and vanishes at loq fi. 

Our results also suggest that the shift in phonon fre- 
quencies has a different dependence on carrier density 
in single layer and many layer systems, eqs. (|24p and 
P5p. For a given carrier density, the shifts in phonon 
frequencies should scale as <Wq L ~ ^Wq^oV^HA-L- As- 
suming that t±/to w 0.1, the shift in a bilayer should be 
smaller than in a single layer sample with the same car- 
rier concentration. This is consistent with experimental 
results which show that the phonon frequencies in single 
layer systems are consistently lower than in samples with 
many layer s 8 ' 9 ' 10 . The difference between a single layer 
and many layer systems is due to the fact that the low en- 
ergy electronic wavefunctions has a reduced weight on the 
sites connected to other layers. Hence, it depends on the 
stacking order, and the shift is different in samples with 
regions with rhombohedral structure, (123123 • • • In 
systems where the charge distribution among the layers 
is not uniform, we expect that the shift of the phonons in 
different layers is also different, leading to a splitting of 
the single layer phonon frequencies. Similar effects may 
also occur in bulk graphite 35 . 

In summary, we have studied the effect of electronic 
inhomogeneities in the phonon spectrum of Raman active 
modes in graphene. We have shown that the electron- 
phonon coupling leads to a shift of the optical phonon 
frequency that is dependent on local electronic density. 
We argue that the frequency shift is larger at the edges 
than in the bulk of graphene and its value is of order of a 
few wavenumbers. These results have their origin on the 
Dirac-like nature of the quasi-particles in these materials 
and hence do not have an analogue in the ordinary 2D 
electron gas. 
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FIG. 2: Notation used for the atomic displacements used in 
the text. 



APPENDIX A: ANALYTICAL MODEL FOR THE 
IN PLANE PHONONS IN GRAPHENE. 

1. The model. 

The simplest model for the phonons of a single 
graphene plane includes only nearest neighbor cen- 
tral forces^, following similar models for the diamond 
lattice^: 



2 

mn 

2M 



E jfofei - a m »)(r fci - r m „) 
2 

k,l\m.n 



( A1 ) 

where the indices fc, I and m, n label lattice sites which 
are nearest neighbors. The three possible orientations 
of the bonds attached to a given site, m, n, allows us 
to define three unit vectors, h l mnl i = 1,2,3. We define 
the displacement of the atom at site m,n, r mn by the 



three projections x\ 



These numbers satisfy 



i 2 3 x mn = 0- The model contains a single parame- 
ter, luq = \J K/M, where K is the spring constant of the 
bonds, and M is the mass of the Carbon atom. 
The equations of motion are: 



U X mn ~ V, 



') 2 ^ ^ { X rrm 



x 3 ,, ,,) 

• h m"n" I 



where the indices 



(A2) 

' and m"n" label sites which are 
nearest neighbors to site ran (see Fig. [5]). We now define 
the variable: 

bmn — x m ' n ' + x m "n" + x m"'ri" (A3) 

using the displacements at the three sites connected to 
site mn (see Fig. [5]). In terms of these variables, the 



equations of motion, eq. (|A2|) , can be written as: 



where 



label the three sites connected to site 



The calculation of the phonon eigenstates is reduced to 
a tight binding model with one orbital per site in the 
honeycomb lattice. From eq. (|A4[) . we obtain two bands: 



1 



■ = J -±-i/3 + 2[cos(k • ai) + cos(k • a 2 ) + cos(k • a 3 )] 
loq V 2 2 

(A5) 

where ai , a 2 are the unit vectors of the honeycomb lat- 
tice, and and a 3 = ai — a 2 . From the knowledge of the 



variables b mn the equations of motion, eq. (|A2[) . can be 



written as: 



,2^.1 



2 [Xk 
3^ 2 



kl 



i ) 2 mn 



(x 1 



(A6) 



From these equations the atomic displacements can be 
deduced from the set {b mn }. 

The equations of motion, eq. (|A2[) assume that all the 
variables b mn are different from zero. When b mn = 0, 
the equations (|A6|) admit two additional solutions for 



uj = 0, and u> — 3wq/2. The bands obtained in eq. (|A5j) 
correspond to the longitudinal acoustical (LA) and lon- 
gitudinal optical modes (LO). The two additional flat 
bands obtained when b mn = describe the transverse 
acoustical (TA) and transverse optical (TO) modes. The 
phonon bands are shown in Fig. [5] . The existence of a flat 
TA band at ui — reflects the band that the honeycomb 
lattice can be distorted without changing the distance be- 
tween nodes. These deformations do not have an energy 
cost in a nearest neighbor central forces model described 
in eq. (|Aip . The velocity of sound of the LO modes is 
v s — (woa)/(2-\/2), where a is the lattice constant. 
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X 2 X 3 




FIG. 4: Displacements used to define the variable QirriTL st an 
atom at a zigzag edge (red), and the variable 6 mn (blue). 



Thus, the equations of motion of the atoms can be 
mapped onto a tight binding model. The only differ- 
ence with the bulk case is that the description of the 
displacements of the atoms at the boundary require the 
definition of two effective orbitals. The position of the ef- 
fective orbital level a m „ at the edge, 0Jq/2, is lower than 
that for the variable b mn , (3wq)/2. This reflects the fact 
that atomic fluctuations are enhanced at the edge. 

Spectra of longitudinal phonons. 256 sites. 



2. Defects 

The mapping to a scalar tight binding model of the 
Hamiltonian in eq. (|Aip can be extended to lattices with 
defects. We describe the defect as the absence of bonds. 
Hence, an atom near a defect is attached to fewer neigh- 
bors than one at the bulk. This implies that the condition 

x mn = is no longer satisfied. We can take this into 
account by defining a new variable at the sites near the 

defect, Cmn — 3^ ram 

bonds which remain intact. 



where the sum is restricted to the 



Phonon density of stales. Vacancy. 
Fluctuations in bond length. Vacancy. 
Phonon density of stales. No vacancy. 
Fluctuations in bond length. No vacancy. 




a. Zigzag 



Frequency 



The atoms at a zigzag edge are connected by only two 
bonds to the rest of the lattice. We define the variable 
b mn using the displacements of the two nearest neighbor 
atoms to the edge atom mn. The atomic displacements 
used to define the variables a mn and b mn are sketched in 
Fig. [3]. The equations of motion for the variables a mn 
and b mn when the indices ran label an atom at the edge 
become: 



2i 
2 i 

U> O m i n i 



2 (^mn bjnji) 
Orn,n - 



2 mn 2 



m'n' 
m"n" 



^a mi (A7) 



where the indices m'n' label the sites which are the near- 
est neighbors of the vacancy, and m"n" stand for the 
next nearest neighbors. The equations of motion for the 
remaining atoms are not changed from eq. (|A4[) . 



FIG. 5: Phonon density of states and spectral function of 
bond fluctuations with and without a vacancy. 



b. Vacancy. 

As in the case of an atom at a zigzag edge, the three 
atoms near a vacancy are connected by bonds to two 
nearest neighbors only. As in the previous case, a new 
variable, a mn = Yl'i x mn needs to be defined at these 
three sites. The equations of motion for the variables 
a mn and b mn are those in eq. (|A7[) . 

The phonon density of states in clusters with and with- 
out vacancies, and the spectral strength of the bond 
length fluctuations, in the bulk and near a vacancy are 
shown in Fig. [5]. Contrary to what happens for the 7r 
electronic band, the phonons are not too disturbed near 
a vacancy, although some shift of spectral strength to 
lower energies takes place. 
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